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1 Introduction 

Let Q C R n be a bounded C 2 domain. Let A : Q — > M n (R) (where M n (M) is the space 
of n x n matrices with real coefficients) be a bounded measurable function satisfying the 
following uniform ellipticity condition: there exists c > such that, for almost every x £ Q 
and all (6l n , 

A(x)Z-£>c\t;\\ (1.1) 

where, for all u,v £ M™, u ■ v is the standard scalar product of u and v and |£| stands for 
the Euclidean norm of £. Given / £ W~ 1,2 (Q), the Lax-Milgram theorem shows that there 
exists a unique function u £ Wq ,2 (VL) such that 

div(AVw) = /, (1.2) 

in the sense that, for all function tp £ T>(Q), 

I A(x)Vu(x) ■ V(p(x)dx = -</, <p). (1.3) 
Jn 

A classical result of N. Meyers (p2j) asserts that there exists e > 0, only depending on Q 
and A, such that, for all p £ (2, 2 + e), if / £ W -1 *^), then it £ iy o 1,p (fi) and 

ll^ll VK 1 ^(Q) — C p ||/|| w-i.P(n) 

where C p > only depends on Q, A and p. This results still holds when Q is Lipschitz (|llj). 
Let us now focus on the case when n = 2 and 11 C M 2 is a polygonal domain. Let A 
as before. We shall be interested in approximate solutions of (11.2j) obtained by Galerkin 
approximations based on spaces of piecewise polynomial functions. Let us briefly present 
these approximation schemes. By triangulation of Q, we mean a family T of triangles such 
that 

n={jT. 

TeT 

Say that T is admissible if and only if, for all T, T' £ T, T D T" is either empty, or a vertex, 
or a common side. If T is an admissible triangulation of Q, for all T £ T, denote by /it 
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the diameter of T and by pr the inner diameter of T {i.e. the diameter of the largest ball 
included in T). Define also 

V\T) := {v G C(fi); u| T ePiVTe T} 

where Pi denotes the space of polynomials with degree less or equal to 1. Define also 
Vq(T) as the subspace of V l (T) made of functions vanishing on dQ. It is plain to see that 
V\T) C W l ' 2 {tt), so that Vq (T) C W^iVt). Another application of the Lax-Milgram 
theorem gives a unique function uj- G Vq(T) such that 

/ A(x)Vu T (x) ■ Vv(a;)(ix = - / f(x)v(x)dx \fv G ^(T) (1.4) 
in Jn 

(note that the bilinear form given by the left-hand side of (11.41) is still coercive on Vq(T) x 
Vq(T)). Assume now that, for all h > 0, we have an admissible triangulation Th of f2, such 
that 

sup fyp = h. 
TeT 

Say that the family of triangulations (7h) h>0 is regular if and only if there exists a > such 
that 

— < cr VT G 7fe, V/i > 0. 
For all /i > 0, write tt^ instead of Uj- h . Then, one always has (see [5], Chapter 5): 

lim|K -u\\ wh2{a) = 0. 

In the present work, we establish a Meyers type result for the approximate solutions Uh of 
(II. ip . as well as an improved convergence result of to u in a Sobolev norm. More precisely, 
we establish: 

Theorem 1.1 Assume that Q C M 2 is a convex polygonal subset and that A G L°°(Q) is 
uniformly elliptic, in the sense that A satisfies (11.11) . Then, there exist e > and C > such 
that, for all p G (2, 2 + e) and all f G W~ l ' p (VL), the solution u h G Vq(%) of (Ol) belongs 
to W 1 *^) and, for all h>0, 

\\ u h\\wi>P(Q) — C ll/llw- 1 -?^) • (1-5) 

Moreover, for all p G (2, 2 + e), 

]im||u fc -u|| wl . p(n) =0. (1.6) 

Theorem 11.11 is optimal in the scale of Sobolev spaces in the following sense: 

Proposition 1.2 Assume that Q C M 2 is a convex polygonal subset. Let e > 0. Then, there 
exists A G L°°(Q) uniformly elliptic with the following property: for all p > 2 + e, there 
exists f G Vy~ 1,p (f2) such that, if Uh G Vo(Th) is the solution of (11.41) . the (uh)h>o are n °t 
uniformly bounded in W 1,P (Q). 



3 



As an immediate consequence of Theorem 11.11 and the usual Sobolev embeddings, we get a 
version of Theorem 11.11 for Holder spaces: 

Corollary 1.3 Under the assumptions and notations of Theorem \l.l\ there exist r] > and 
C > such that, for all pe(2,2 + e) and all f G W'^iVL), 

\\ u h\\co<n(ji) — C WfWw-^pin) (1-7) 

and 

lim \\u h - u\\ c0 , v(n) = 0. (1.8) 

When / G L 2 (f2) and if one furthermore knows that the solution u of (11. II) belongs to 
VF 2 ' 2 (f2) (this is the case for instance if A G C 1 (f2)), one proves the following conclusion: 

Corollary 1.4 T/ie assumptions and notations are the same as in Theorem Assume 
that, for all f G L 2 {Vl), the solution u of (11.11) belongs to W 2 ' 2 (£l). Then, for allp G (2, 2+e), 
there exist 9 G (0, 1) and C > such that 

IK - Hlun^) < ll/lli^n) • (1-9) 
Let us give a few comments about these results. 

First, at fixed h, Uh is piecewise linear. Moreover, since Uh is actually given by a finite system 
of linear algebraic equations, it is easily seen that there exists Ch > such that 

IKII w i,<» ( n) ^ °h \\f\\i?(n) ■ ( 1 - 10 ) 

However, as Proposition II. 21 shows, the constant Ch in ( 11.101) cannot be bounded when h — > 
in general. 

Theorem 11.11 and Corollary 11.31 can be seen as regularity results for approximate solutions of 
(11. 2p obtained by the finite elements method. These results are the counterparts of classical 
Meyers and De Giorgi type regularity results for the solution of dw(AVu) = f (see the 
results in [TT], [T7j stated at the very beginning of the paper). Since Q C 1R 2 , Theorem 11.11 
yields Corollary 11.31 at once by means of the classical Sobolev embeddings. 
Under the extra assumption that A is symmetric and when the data / belongs tp L 2 (Q), a 
version of Corollary 11.31 was proved in [19], Theoreme 6.2. Here, we get rid of the symmetry 
assumption on A and we also derive the convergence result in (jl.8p . while only (II. 7p was 
established in [Tj5]. To our best knowledge, Theorem 11.11 and Proposition 11.21 are completely 
new. 

Using ideas analogous to those involved in the proof of Theorem 11.11 we also establish 
Sobolev, L°° and C°' v estimates for second order uniformly elliptic operators on quite general 
graphs. A lot of work was devoted to estimates for reversible Markov chains on infinite graphs 
(see for instance [71 El ED] and the references therein). In these references, the operators under 
consideration have positive coefficients. In the present work, we deal with more general 
operators, for which no positivity assumption is made on the coefficients but which are still 
elliptic (see Section [5] below for the statements and the proofs). It seems therefore impossible 
to use techniques such as the Moser iteration or the maximum principle. 
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Let us describe our strategy for the proof of Theorem 11.11 and its corollaries. As in [19] , the 
proofs go through the analysis of some second order elliptic operators on graphs. However, we 
simplify the arguments in a significant way. Given a triangulation Th, we define a graph r\, 
the vertices (resp. the edges) of which are the vertices (resp. the edges) of Th, and consider 
a suitable "second order elliptic" operator on Th, denoted by which approximates L on 
T h . We then reduce (jl.5p to proving that the inverse of L h is L 2 (T h ) — W 1,2+£ (T h ) bounded 
for some e > 0. 

The proof of this boundedness result for L^ 1 is inspired by an analogous perturbation argu- 
ment for 2m- th order operators in IR 2m (see [T], Chapter 1). Namely, we first observe that 
is bounded from W 1,p {Th) to W~ 1,p {Th) for all p G (1, +oo), and that it is an isomorphism 
from W^iTh) into W^ 2 {V h ). Then, using the fact that the W l > p (T h ) and the W~ ljp {V h ) 
form an interpolation scale for the complex method (this result is also new and has its own 
interest), we use a general perturbation result to conclude that is an isomorphism from 
W 1,p (Th) onto W~ 1,p (Th) for p e (2 — e, 2 + e) for some e > 0, which completes the proof. 
We argue similarly to obtain W 1,p , L°° and C v estimates for second order elliptic operators 
on graphs T, under some geometric assumptions on T. In particular, the volume of a ball 
with radius r > has to be controlled from below by r 2 . Sobolev estimates are obtained 
by a perturbation argument relying on interpolation results for Sobolev spaces. Then, the 
geometric assumptions on T imply that L 2 (T) embeds into W~ 1,P (T) and W 1,P {T) embeds 
into C" ? (r) for p > 2, which yields Holder estimates. 

Theorems 11.11 and its corollaries are stated for open polygonal subsets of IR 2 . It is plausible 
that an adaptation of our method would yield an analogue of Theorem 11.11 for polyhedral 
subsets of IR 3 . However, this would not yield an analogue of Corollary 11.31 Obtaining a De 
Giorgi type regularity result for approximate solutions of ( 11. 2ft by a finite elements method 
is an open challenging problem. 

The plan of the paper is as follows. In section [2, we present the discrete setting. Section [3] 
is devoted to general interpolation results for Sobolev spaces on graphs for the real and the 
complex method. In Section HI we show Theorem 11.11 and its corollaries, while estimates for 
general second order elliptic operators on graphs are stated and established in Section [5j 
Notation: if A(f) and B(f) are two quantities depending on a function / ranging in a set 
E, say that A(f) < B(f) if and only if there exists C > such that, for all / e E, 

A(f) < CB(f), 

and that A(f) ~ B(f) if and only if A(f) < B(f) and B(f) < A(f). 
Acknowledgements: we would like to thank Franck Boyer for helpful discussion. 

2 Presentation of the discrete framework 

Let us give precise definitions of our framework. The following presentation is partly bor- 
rowed from [10J. In the sequel, T stands for a non-empty finite or infinite set. When T is 
finite, fix a non-empty strict subset of T, called the boundary of T and denoted by dT. 
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2.1 The metric on T 

Let h xy = h yx > a symmetric weight on Y x T. We assume that h xx = for all x G Y. If 
x, y G T, say that x ~ y if and only if h xy > or x = y . We define for all x G T 



xy 



h x := sup h 
and 

ft, := sup h X y 

x,y£F 

Denote by E the set of edges in T, i.e. 

e ■= {(x, y ) e r x r ; x~ y }, 

and notice that, due to the symmetry of h, (x, y) G E if and only if (y, x) G E. In the sequel, 
we always assume that V is locally uniformly finite, which means that there exists iV G N* 
such that 

sup e r ; y ~ x} < N (2.11) 

(here and after, §A denotes the cardinal of any subset A of Y). 

We also assume that h is locally uniformly controlled with a constant C\y > 1, which means 
that for all x, y, z G Y such that x ~ y with x ^ y and x ~ z with x ^ z, we have 

<^<CV- (2.12) 



CV ft__ 

For x, y G T, a path C joining x to y is a finite sequence of vertices x = x, x N = y such 
that, for all < i < iV — 1, Xj ~ Xj + i. By definition, the length of such a path is 

JV-l 



k=0 

Assume that Y is connected, which means that, for all x, y G Y, there exists a path joining 
x to y. For all x, y G Y, the distance between x and y, denoted by d(x,y), is defined as 
the infimum of the lengths of all paths joining x and y. For all x G T and all r > 0, let 
B(x,r) := {y G Y, d(y,x) < r}. 

Remark 2.1 If x G T and r > are st/cn taat r < taen B(x,r) = {x}. Indeed, let 
y G -B(x, r) ; so tnat d(y, x) < r < If y ^ x, there exists z ~ x snca i/iai 

/ice 

^ < d(y, x) < — < h xz , 

which is impossible. 

If B = B(x, r) is a ball, set aB := B(x, ar) for all a > 0. 
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2.2 A measure on the graph and the edges 

We also consider measures on E and T. Assume that, for all x ~ y G T, we are given a 
coefficient \i xy > such that fi xy = n yx . Assume also that there exists C > such that, for 
all x, y, z G T satisfying x ~ y and x ~ 2, 



For all a; G r, set 



^ < C. (2.13) 



m m 



if Act, let 



m(A) := 2_^m(x). 



For all 1 < p < +00, say that a function / on T belongs to L P (T, m) (or L P (T)) if and only if 

\ Vp 

^|/(x)| p m(x) <+oo. 
.xer / 

Say that / G L°°(T, m) (or L°°(r)) if and only if 

ll/lloo :=sup|/(x)| <+oo. 

When T is finite and 1 < p < +00, say that / G Lq(T) if and only if / G L p (r) and f(x) = 
for all x G <9r. 

Remark 2.2 It zs an elementary observation that the space of functions on T with finite 
support is dense in L P (T, m) for all 1 < p < +00. 

For all 1 < p < +00, say that a function F on E belongs to L P (E, fi) (or L P (E)) if and only 
if F is antisymmetric (which means that F(x,y) = —F(y,x) for all (x,y) G E) and 

ll*V= ( E \F(x,y)\ P ^y) <+oo. 

\(*,2/)e£ / 

Say that F G L°°(E, fj,) (or L°°(E)) if and only F is antisymmetric and 

ll^lloo := SU P l-^fosOl < +°°- 
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2.3 Differential and gradient 

In the spirit of [3], for all function / : T — > R and all (x, y) G E, define 

df(x, y) = I K, ltxfJ - 
y if x = y. 

Say that d is the differential operator on Y. 

We also define the sublinear operator "length of the gradient" by 

1/2 



v ^) = /rfe l/(?/) " /(2:)|2 ) 

x \y~x J 



for all function / on T and all x G T. It is a consequence of (12. lip . ( I2.12p . (I2.13P and (I2.14p 
that, for all function / : T — > R and all p e [1, +oo], 

\\df\\LnE,»)~\\Vf\\L H r im )- 
2.4 Sobolev and Holder spaces 

Let 1 < p < +oo. Say that a scalar-valued function / on T belongs to the (inhomogeneous) 
Sobolev space W^iY) (see also [31 M, W if and only if 

ll/ll W i.p(r) := \\f\\ L p(r) + ll v /llLP(r) < +°°- 

In the case when T is finite, denote by W^'^tT) the subspace of W 1,P (T) made of functions / 
such that / = on dT. 

It is then routine to check that both W^ P {T) and W 1,p (r) are Banach spaces. 

Finally, if 1 < p < +oo, as in the Euclidean case, we define W~ 1,P (T) as the dual space 
of W 1,p '(r) (when T is infinite) or of Wq' p (T) (when T is finite), equipped with its natural 
norm. Here and after, p' is defined by - + -7 = 1. 

Remark 2.3 It is easy to check that, for all 1 < p < +00, the space of functions with finite 
support is dense in W 1,P (T). 

For all fi > 0, we define the discrete homogeneous (resp. inhomogeneous) Holder spaces of 
exponent fi as 

:= {u E R r ; M CM(r) = sup K * ) ~^ )l < +00} 

and 



x ^y d(x,yY 



C M (r) := {u G M r ; ||w||o(r) = ll M IU°°(r) + \ u \c»(r) < +°°}- 
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2.5 Geometric assumptions on T 

Let us now introduce some geometric assumptions on (T,m). 

1. Say that (T,m) satisfies the doubling property if there exists C > such that, for all 
x G r and all r > 0, 

V(x,2r) < CV(x,r). (D) 

Note that this assumption implies that there exist C, D > such that, for all x E T, 
all r > and all 9 > 1, 

V{x,6r) < C9 D V(x,r). (2.15) 



Remark 2.4 Observe also that, when T is infinite, it is also unbounded (since it is 
locally uniformly finite) so that, if l[D\) holds, then m(T) = +oo (see [16]). 

2. Say that (T,m) satisfies (L a ) if there exist a > 0, c > such that, for all x G T and 
all r > 0, 

V(x,r)>cr a . {L a ) 

3. If 1 < q < +oo, say that (T,m) satisfies a scaled L q Poincare inequality on balls (this 
inequality will be denoted by (P q ) in the sequel) if there exists C > such that, for 
any x G T, any r > and any function / on T, 

E \f(y)-fB\ 9 ™(y)<Cr« lv/(y)IMy), (n) 

y£B(x,r) y£B(x,r) 

where ^ 

/b = TTTm f( x ) m ( x ) 

^ ' x£B 

is the mean value of / on B. 

When T is finite, we will consider local versions of these properties. Let r > 0. Say that V 
satisfies (Di oc ) if there exists C > such that, for all < r < ro and all x G T, 

V(x,2r) < CV{x,r). 

Say that Y satisfies (L a ^ oc ) if there exists C > such that, for all x G T and < r < r , 

y(x,r) > Cr a . 

Let 1 < g < +oo. Say that T satisfies (Pqjoc) if there exists C > such that, for all x G T, 
all < r < r and any function / on T, 

\f(y)-fB\ q m(y)<Cr« \Vf(y)\ q m(y). 

y&B(x,r) y£B(x,r) 
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3 Interpolation and embeddings of Sobolev spaces 

3.1 Real and complex interpolation of Sobolev spaces: the case of 
unbounded graphs 

Let us first focus on the case when T is infinite, and therefore unbounded. We will make 
constant use of the following notation here and in the sequel: 

Definition 3.1 When (P q ) holds for some 1 < q < oo, we defineqo = inf {q G [1, oo) : (P q ) holds}. 

Note that if (P q ) holds for some q > 1, then q < q (see [15]), and that (P r ) holds for all 
q < r < +oo. 

Let us first state the following real interpolation theorem for Sobolev spaces: 

Theorem 3.2 Let q G [1, +oo) and assume that ID]) and (P q ) hold. Then, for all 1 < r < 
q < p < +oo, W^iY) = (W l > r (T), W 1 ' 00 ^))^ _. 

The corresponding result for Sobolev spaces on Riemannian manifolds was proved in [2], 
Theorem 1.1. We established in [3], Theorem 1.18, a statement very similar to Theorem 13.21 
for homogeneous Sobolev spaces, and the proof, which we skip here, is analogous to the one 
of Theorem 1.18 in [3], even if the definition of the gradient is slightly different. The only 
important assumptions to be able to interpolate are (ID]) and (P q ), see Theorem 7.11 in [2]. 

As an immediate corollary, we obtain: 

Corollary 3.3 (The reiteration theorem) Assume that T satisfies ID]) . (P q ) for some 
1 < q < +oo. Let q < pi < p < p 2 < +oo, with - = 1 , then 

P Pi P2 

w lj, (r) = (w^n.w'^cr)) . 

If qo — 1, then one can also take 1 = p\ < p < p 2 < +oo. 

Let us now turn to the complex interpolation of Sobolev spaces. We first recall the following 
complex reiteration theorem: 

Theorem 3.4 \9j/ For any compatible couple of Banach spaces (Ai,A 2 ) we have 

(A 1 ,A 2 ) PtP 

for all Ai, A2 and a in (0, 1) and all p\.p 2 in [1, +00], except for the case p\ = p 2 = 00. Here 
(3 and p are given by (3 = (1 — a)Ai + a\ 2 and i = + ^. 

From this, we deduce a complex interpolation result for Sobolev spaces: 

Corollary 3.5 Assume that T satisfies ID]) and (P q ) for some 1 < q < +00. Let qo < P\ < 

p < p 2 < +00 where qo is given by Definition ^. 1\ Then for a = V__X_ = p\^ 2 Zpl] one has 

[W 1 ' pl (T),W 1 ' P2 (T)] a = W^ P {T). 
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Proof: We apply Theorem GLl with A x = W 1 ' 1 ^), A 2 = W l '°°(T), A x = 1 - i, A 2 = 1 - i 
and 13 = 1 — - . We obtain 

[^(r), w^(r)] Q = (^(r),^ 1 ' 00 ^))^ 

and the result follows from Theorem 13.21 which yields 

since p > q . □ 
Now for dual spaces, we recall the following theorem: 

Theorem 3.6 (UN, Theorem 3.7.1) Let Ai,A 2 a compatible couple of Banach spaces. Then 
for all 9 G (0, 1) and 1 < s < oo, we have 

\(A u A 2 y s y = (Al,A* 2 ) e , s , 
Combining this Theorem with Theorem 13.2} we deduce: 

Theorem 3.7 Assume that V satisfies [D]) and (P q ) for some 1 < q < +oo. Then, for all 
l<ri<r<r 2 <oo with r' > qo, we have 

(w- l > ri {r),w- l > r2 {r)) 6jr = w- l > r (r) 

where 9 G (0, 1) such that - = — + — . 

Proof: We have 1 < r 2 < r' < r[ < oo with r' > q and -V = 4- + ^—r-. Then 

r r 2 r 1 

□ 

Theorem 3.8 Assume that Y satisfies (TP]). (P g ) /or some 1 < q < +oo. Let 1 < pi < p < 
p 2 < q' Q , where qo is given by Definition ^. 1\ Then 

[W- 1 ' Pl (T),W- 1 ^(T)] a = W-^ P (T) 
i i 

where (5 = 1 — - and a = ■ 

PI P2 

Proof: We have p\ < p < p 2 < q' Q . Then Theorem 13.71 shows that 

Pi ,y P2 

and 

(w- 1 ' 1 (r),w- 1 ' oo (r)) 1 _i_ = w- l ' p {v). 

This, with Theorem 13.41 yields 

[W- 1 ' P1 {T),W- I ' p2 {r)] a = (W- 1 ' 1 ^),]^- 1 ' 00 ^))^ = w-^iv). 
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3.2 Real and complex interpolation of Sobolev spaces: the case of 
bounded graphs 

Assume now that the graph T is bounded. 

Theorem 3.9 Let q e [l,+oo) and assume that (Di oc ) and (Pqioc) hold. Then, for all 
l<r<q<p<+oo, W^(T) = {W^^W^))^ . 

Proof We refer to [2] and [5] for the proof. We just mention the differences between the 
proof in [21 |3] and the present situation. We have to estimate the functional K of real 
interpolation. Let 

n = {xeT;M((\Vf\ + \f\y)(x)>a«(t)} 
where M. is the uncentered Hardy-Littlewood maximal function, that is 

Mh(x) := sup — -— V \h(y)\m(y) 

B3 X V(B) ^ 

and a(t) = M((\Vf\ + |/|) 9 )*i(t).We recall that for a function g on T, g* denotes its 
decreasing rearrangement function. 

If fl = T (note that this may happen since m(T) < +oo), this is the easy case, we argue as 
in |2], Section 3.2.2. 

Otherwise, we write T as the union of balls B l , of radius p small enough (namely < p < r ) 
and having the bounded overlap property. We take (tp l )i a partition of unity subordinated 
to the covering (B' l ) i and write / = ^ f(p l = f l . We estimate the functional K as in 
[2], Section 3.2.2, using the Calderon-Zygmund decomposition applied to /\ Note that here 
the functions (&*)_. of the decomposition should belong to ^^^(r), that is bj = on dT. For 

this, we take b % - = f % x) if B l - fl dT ^ and, if not, we define as usual &*■ = (f l — (Db^X)- 
We characterize then K as in [2]. Integrating K, we get the interpolation result. □ 



As in section 13.11 from this theorem, we deduce: 



Corollary 3.10 Assume that T satisfies (Di oc ) and (P q .ioc) for some 1 < q < +oo. Let 

_i i 

q < pi < p < p 2 < +oo where q is defined by (13. ip . Then for a = = p ^ ^ one has 

[W^(T),W^(T)] a = W^(T). 
As far as dual spaces are concerned, we have: 

Theorem 3.11 Assume that T satisfies (Di oc ) and (P qt i oc ) for some 1 < q < +oo. Then, 
for all 1 < ri < r < r 2 < oo with r' > q , we have 

where 9 E (0, 1) such that - = — + — . 
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Finally, 

Theorem 3.12 Assume that T satisfies (D loc ), (P g j oc ) for some 1 < q < +00. Let 1 < p 1 < 
p < P2 < q'o- Then 

[W- 1 > pl (T),W- 1 > P2 (T)] a = W^iY) 
j 1 

where (3 = 1 — - and a = j 1 r 

3.3 Sobolev embeddings 

Assuming (TDl) . (P g ) and (L CT ) for some cr > 0, we obtain Sobolev embeddings analogous to 
the Euclidean situation. Namely: 

Proposition 3.13 Assume that F satisfies ID]) . (P q ) and (L a ) for some a > 0. Then 

ll/||jy*(r) < ||V/||LP(r) 

for all qo < p < a where p* = -^-^ . 

This result is proved in [20] , Corollary 3.3.3, p. 75. 

Proposition 3.14 Assume that T satisfies (fDj). (P q ) and (L a ) for some a > 0. Then, for 
all p > max(go, o~), 

W X ' P {Y) ^ c"(r) 

with rj := 1 — ^, where the symbol °->- means that the embedding is continuous. 

Proof: Let u e W 1,P (T), x,y G T with x ^ y and set r := d(x,y). Define := B(x,r), 
so that 5 C B(x, 2r) H 2r). One has 

\u(x) - u(y)\ < \u(x) - u B \ + \u(y) - u B \ , (3.16) 

and we estimate each term of the right-hand side of (I3.16p . 
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Let us first focus on the second term. Using Remark 12.11 write 
o 

\ u (y) — u b\ < ^2 \ U B(y,2 k r) ~ u B{y,2 k + 1 r)\ + \uB(y,2r) ~ Ub\ 
k=—oo 

1 1 

k=-oo ^' > zGB{y,2 k r) ^ ' z&B 

° y(y 2 k+1 r) 1 
^ E V(tW) V(y,2^r) 22 \<z)-u B{y ^ r) \m{z) 

k=-oo Vi " ; Vi " ' zeB(y,2 k +^r) 

c 

+ y, 2r ) l«(*)-«B( Wl 2r)K*) 

E IW™w) + Cr(y^ E lv«W" 

fc=-oo y Vy ' ; zeB^^fe+ir) / \ W ' ' z<=B{y,2r) 

< C ( ] || V« || LP( r) + CV 1 "? ||W|| iP(r) 



\k=— oo 



< CV 1 p || Vw || L p ( r) = Cd{x,yf * ||Vw|| ZP(r) , 

where we used (D) and (P p ) in the third and fourth inequalities and (L a ) in the fifth one. 
Arguing similarly, one obtains 



\u\x\ 



u B \ <Cd{x,y) 1 ^ \\Vu\\ LP[T) . (3.17) 



Thus u G ^(r) and Hc^p) < C\\ Vu \\lp(f)- 

It remains to prove that u G L°°(T). For all x G T, write 



(x)| < |w(x) - us^.i)! + \u B (x,i) - u(y) \ + \u(y)\ (3.18) 



for all y G B(x,l). 

Applying ( 13.1 7p with r = 1, one obtains 

- U B (x,l)\ < C \\^ U \\lp(B(x,1)) ■ ( 3 - 19 ) 

Then, taking the average over B(x,l) in (I3.18P yields 

\u(x)\ < \u(x) - U B (x,l)\ + y} E -MB(x,l)[^(y) + y / ^ l U (^)l m (3/) 

~ ° I 1 + t/v ^ ) 11 Wu Wlp{b{x,x)) + — — -rlMkcBM)) 

\ K (X, 1)P J v[X, 1)p 

< C\\u\\ w i, P (r). 

We used (13.191) . (P p ) and (L a ) in the last inequality. Thus ||M||L°°(r) < CIMI w^r) and 
therefore u G (7^(17) with ||u||c(r) < C|MI vT^r)- □ 

14 



4 Meyers type regularity for Galerkin schemes 



This section is devoted to the proof of Theorem 11.11 Proposition II. 2\ Corollary 11.31 and 
Corollary Ol 

4.1 Definition of a graph and an operator 

Let Q C 1R 2 be a convex polygonal domain and {Th) h>0 be a regular family of triangulations 
of Q as in the statement of Theorem II. 11 Let T/j be the set of vertices of Th- Define 
dT h := T h fl 80. . For all x, y G T h , say that x ~ y if and only if x = y or there exists a 
triangle T in Th such that the segment joining x and y is an edge of T. For all x ~ y G T^, 

set 

where we recall that \u\ denotes the Euclidean norm of u G JR 2 . Note that, since the 
triangulation is admissible, (12.111) and (12 . 1 2p hold. For all x ~ y, define also 

so that (12.131) holds and 

m(x) ~ h\ 

for all x G Th- 

Lemma 4.1 The graph Th, endowed with h xy and \x xy satisfies (Di oc ), (P2,ioc) and (L2,/ c); 
with constants only depeding on Q and a and independent of h. 

Proof: The proof uses essentially the regularity of the graph Th and Poincare inequality 
{P%ioc) on ^- A detailed proof can be found in [19], Chapter 7, section 7.2. □ 

For all i e Tj, fl fl, define (p x as the unique function in Vo(7h) suc h that 

for all y G Th- The functions (<p x )xeT h nn form a basis of V^Th)- 
Let Rh be the "reconstruction" operator defined by 

RhUh ■= ^2 
xer h nn 

for all function Uh on T h vanishing on dT h - In other words, RhMh is the function in V^(7/i) 
which coincides with Uh at the vertices of Th- We claim that the norms of Uh and RhUh in 
the various functional spaces introduced before are equivalent: 

Proposition 4.2 Let p G (1, +oo) and r\ G (0, 1). For all function Uh on Th vanishing on 
dT h , 
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2 - \\Vuh\\ L v(r h ) ~ \\V R hUh\\ LP{n) , 

3. ||Wft|lc"7(r h ) ~ \\ R hUh\\co,v(u)- 
The implicit constants in these equivalences of norms only depend on a and are independent 
ofh. 

Proof: It is a straightforward consequence of the definitions of the norms, of i? ft , the 
properties of the triangulation: namely the regularity of T ft and the fact that it is locally 
unifomly bounded, and finally the Poincare inequality on a bounded open convex subset of 
IR 2 . A detailed proof can be found in Appendix E of |19j . □ 
Let us now define a maximal accretive operator on L^{Th) associated to the Galerkin scheme. 
For all Uh,Vh G WQ ,2 (T h ), define 

Qh {uh, Vh) ■= / A(x)V {R h Uh) (x) ■ V {R h v h ) (x)dx. 
Jn 

The operator Vi? ft : Ll(T h ) — > L 2 (Q,C n ) is defined on the whole L^(Y h ) space, and the 
operator F \-t AF is obviously bounded on L 2 (fl,M. n ). Moreover, if u~h G W 1,2 (r ft ), 

A(x)V {Rku~ h ) (x) ■ V {RhUh) (x)dx > 8 || V (R h uQ \\ 2 L2{n) 

> \\R h u h \\ 2 L2{n) 

As a consequence, using a construction due to Kato (p3]), there exists a unique maximal 
accretive operator L ft on L^iTh) such that, for all Uh G T>(Lh) (here and after, T>(L) stands 
for the domain of an operator L) and all Vh G W ,2 (Th), 



(L h u h , v h ) L 2, Vh) = / A(x)V (R h u h ) (x) ■ V (R h v h ) (x)dx. 
Jn 

The following facts are easily checked: 

Proposition 4.3 For all p G (l,+oo) ; the operator Lh extends to a bounded operator from 
Wa' p (r h ) to W'^iTh). Moreover, Lh is an isomorphism from W ' (r ft ) into W 1,2 (T ft ). 
The norms of Lh and of its inverse (forp = 2) do not depend on h. 

Proof: Let p G (1, +oo) and p' G (1, +oo) such that ± + ^ = 1. Let u~h G W 1,p (r fe ). The 
map v\ i->- J n A(x)V (RhUh) (x) ■ V (RhVh) (x)dx is clearly a linear functional on Wq' p (r h ). 
Moreover, for all Vh G Wq' p (r ft ), 

\(L h Uh,v h )\ < \\V(R h u^)\\ LP{n) jV(R h Vh)\\ LP ' (n) 

This ends the proof of the first part of the claim. Assume now that p = 2. For all 
bounded linear functional T on W ' 2 (Th), the Lax-Milgram theorem yields a unique func- 
tion G W 1,2 (r ft ) such that, for all v h G W 1,2 (r h ), (L h Uh, Vh) L 2 {r h ) = T(v h ), and one 
has IIm^IIvjai. 2 ^) ^ H^lliy- 1 . 2 ^)- This shows that L h is an isomorphism from W^iTh) into 
W~ l,2 (Th)- Moreover, the norms of L ft and L^ 1 do not depend on h because of Lemma [4.11 
and Proposition 14.21 □ 
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4.2 Conclusion of the proof 



Our next step is a perturbation argument, in order to show that is still an isomorphism 
from WQ ,p (Th) into W~ 1,p (Th) for p close enough to 2. Indeed, since (P2,ioc) holds, there 
exists qo < 2 such that, for q < p < p < p\ < q' , the W 1,P (T) and W~ 1,P (T) form a 
complex interpolation scale (see Corollary 13.51 and Theorem 13. 8j) . We then combine this fact, 
Proposition 14.31 and the following general result (|21j): 

Lemma 4.4 Let X s , Y s , s G [0, 1] be two scales of complex interpolation Banach spaces. If 
T : X s — > Y s is bounded for each s G [0, 1], then the set of s for which there exists C > 
such that \\Tf\\Y» > C\\f\\x" holds for all f G X s is open. 

Thus, there exists e > such that Lh is an isomorphism from Wl' p {Th) into W~ l ' p iTh) for 
all p G (2 — e, 2 + e). Thus, we have established: 

Proposition 4.5 The operator L^ 1 maps W~ 1,p (Th) into Wl' p (Yh) forp G (2,2 + e) with a 
norm independent of h. 

Let us now conclude the proof of Theorem ll.il Take / G W~ 1,P (VL) and define, for all x G Th, 




We claim that, if Uh G ^(T^) solves (II. 4p and Uh G WQ ,2 (T h ) is defined by 




when i G Td H SI, 
when x G Th D 




(4.20) 



Indeed, L h Uh G 1,2 (r/ l ), and f)4.20p exactly means that, for all x G T h fl Q, 



(L h u h ,ip x ) = (fh,Px) 



where Lp x is the function in WQ 2 (T h ) defined by (p x (y) = 5 xy for all y G T h , so that <p x = Rh^. 




and 




(4.22) 



(/, <Px) 

m{x)Jh{x) 

(fh, fx), 
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where we used (11.4j) in the third equality and the definition of fh in the fourth one. Gathering 
(IQTfl and (jHHD , we obtain ( IQUj) . It therefore follows from ( H~2Uj) and Proposition SSI that, 
for all pe (2,2 + e), 

ll u 'i|lw :l 'i'(r h ) ~ llAllw" 1 .?^) • 
Proposition 14.21 shows that, for all p G (2, 2 + e), 



Moreover, 

II All w-i,j>(r h ) — C 11/11 vK- 1 .p(n) 
Indeed, let f be a function in with finite support. Then, 



(4.23) 

(4.24) 
(4.25) 



\(fh,v)\ 



< 



^2(f, tp x )v(x) 

x&r h 
x&r h 



< 



x& h 

lw/-i,p(c) ll-^^llw^cn) 
which proves ( f£25jl . Gathering (T4T231 . KMf and fT4T25|) yields (fl~5l) . 

For the proof of ( II .6|) . it is enough to interpolate between ( II .5|) and the fact that 
ll w ~~ u ^llwi.2(n) when /i — > 0. □ 
Let us now prove Proposition 11.21 Let e > 0. By [17], there exists a uniformly elliptic 
matrix A G L°°(fi) with the following property: for all p > 2 + e, there exists / G W~ 1,P (Q) 
such that the solution of — div(AVw) = / in Wq 2 (Q) does not belong to W 1,P (Q). Assume 
now that the Uh are uniformly bounded in W 1,P (Q). Up to a subsequence, Vuh converges 
weakly in L P (Q) to some vector field V. But since Vuh —> Vw in the L 2 norm, it follows 
that Vw = V and that V« G L p (f2), a contradiction. □ 
As explained in the introduction, Corollary 11.31 follows at once from Theorem II. II by Sobolev 
embeddings. Finally, for Corollary 11.41 it is well-known (see [5], Chapter 5) that, in this 
situation, 

\\u - u h \\ wl , 2{n) < Ch ||/|| x2(n) . (4.26) 
By the results of [11] stated in the introduction, and since L 2 (Q) <^-> W~ 1,2+e (Vt), one has 

\\u\\ wl ,, +Hsl) < C \\f\\ Lm . (4.27) 

so that 

\\u - u h \\ wl , 2+s{n) < C ||/|| £a(n) • (4.28) 
Interpolating between (14.281) and (I4.26P gives (11.91) with 



1 

2+E 

1 

2+e 
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5 Estimates for general second elliptic operators on 
graphs 

In the present section, we prove estimates in L°° and Holder spaces for general second order 
elliptic on graphs, assuming suitable geometric properties on the graph T, which is assumed 
to be unbounded. 

5.1 Definition of the operators 

Assume that, for all x ~ y, a coefficient c xy G C is given, and that there exist > such 
that 

\c xy \ < Coo Vx ~ y (5.29) 
and 5 > such that, for all u G W 1,2 (V), 

Re ^c xy \du(x,y)\ 2 /j, X y > 5^2\du(x,y)\ 2 fi xy . (5.30) 

x~y x~y 

We associate to these coefficients an operator, which is the discrete version of second order 
uniformly elliptic operators in divergence form in W 1 . To that purpose, we use again a 
classical construction of maximal accretive operators due to Kato Q14J). The operator d : 
L 2 (r) — y L 2 (E) is densely defined (see Remark 1 2. 2 j) and closed, and the operator F y cF is 
clearly L 2 (E) bounded by (|5.29p . Since (I5.30p holds, there exists a unique maximal accretive 
operator L on L 2 (T) such that, for all u G V(L) and all v G W l > 2 (T), 

Lu(x)v (x)m(x) — c xy du(x,y)dv(x,y)fi xy (5.31) 

xeT (x,y)£E 

(recall that T>(L) stands for the domain of L). 

Remark 5.1 Note that assumption ( 15.301) is satisfied in particular when the coefficients 
c xy are real-valued and satisfy c xy > 5 whenever x ^ y. However, we do not make this 
assumption in the sequel, and all that is needed is (I5.30p . 

We say that L is a second order uniformly elliptic operator on L 2 (T) with ellipticity constants 
Cqo and 5. Recall that T>(L) is dense in L 2 (T) and that L generates a holomorphic semigroup 
on L 2 (T). Note finally that, for all u G V(L), 

Re (Lu,u) > c5\\Vu\\l. (5.32) 

Remark 5.2 Let I < p < +oo. Assume thatue V(L)f\W l > p (T) andv G W 1 ' 2 (T)f]W 1 ' p ' (T) . 
From the definition of L, we deduce that 

\(Lu,v)\ < C||VM|| LP (r)||Vv|| LP ' (r) . 

This and Remark \2. 3\ show that Lu extends to a bounded linear form on W l,p ' (T), with norm 
less or equal to C||u||iyi.p(r) ■ Thus L is a bounded operator from W 1,P (T) to W~ 1,P {T), with 
bound less or equal to C . 
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5.2 The estimates on L 

Our first estimates deal with the resolvent of L. Let u E (0, ~) be such that L is a;-accretive 
on L 2 (T), pi E (f , 7r — u) and 

:= {2 G C; 2 = or |arg z\ < fi} . 
Let / E L 2 (T) and A E £ M . We are interested in weak solutions of the equation 

Lu + \u = f. (5.33) 
By "weak solution", we mean a function u E W 1,2 {T) such that, for all v E W l,2 (T), 

c xy du(x, y)dv(x, y)fi xy + A u{x)v{x)m{x) = f {x)v{x)m{x) . 
(x,y)eE xer xer 

We establish that, if (D), (L 2 ) and (P 2 ) hold, given / E L 2 (T) and A E fQ3l has a 
unique weak solution in H /1 ' 2 (r), and that this solution actually belongs to C V (T) for some 
rj > 0: 

Theorem 5.3 Le£ T be a graph satisfying (D), (L2) and (P)- Let L be a second order 
uniformly elliptic operator on L 2 (T) with ellipticity constants and 5. Let u be such that 
L is u-accretive on L 2 (T) and /i E (|,7r — tu). Then, for all f E L 2 (T) and all A E S M , 
there exists a unique solution u E W /1,2 (r) of H 5 . 3 3 j) . Moreover, there exist C, r/ > only 
depending on the constants in (D), (L2) and (P2) and the constants of ellipticity of L such 
that 



|M| L oo (r) < C |A| - \\J || L 2 (r) 

and 

H cv{r) <c\\\i\\\-^\\f\\ L2{r) . 

Theorem 15.31 is a version of the De Giorgi regularity theorem for elliptic equations in this 
context. 

We also obtain estimates on the kernel of the semigroup generated by L. Before stating the 
Theorem, we introduce the following notations: for all x, y E T, h*_+ y is equal to if x = y 
and is the supremum of the weights of edges where at least one of the two vertices belongs 
to the ball B(x, d(x, y)), when x ^ y. We then define 

Theorem 5.4 Let T be a graph satisfying (D), (L2) and (P)- Let L be a second order 
uniformly elliptic operator on L 2 (T) with ellipticity constants and 5. We then have: 

1. The operator L generates a holomorphic semi-group {e~ tL )t>o on L 2 {T) which has a 
kernel, denoted by K t (x, y) in the sequel. For all u E L 2 (T) and all x ET, 

(e~ tL u)(x) = S £ j K t (x,y)u(y). 

Moreover, there exist C, C, f3 > depending only of the geometric constants and that 
of ellipticity, such that for every (x,y) E T 2 and t > we have 
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a. ift<C'h* xy d(x,y), 

b. ift>C'h* xy d(x,y), 



\K t (x,y)\ < —e ~ T ^~; 



\K t (x,y)\<^e-^. 



2. There exist constants C" and rj > depending only on the geometric constants and 
that of ellipticity such that for every t > and for every x, x', y G T, we have 



\K t (x,y)-K t (x',y)\<^( 



C" fd(x,x') 



/\\V 



5.3 Proofs 

We first focus on Theorem 15.31 We begin with the following observation: 

Proposition 5.5 Let T as above and L a second order elliptic operator on T, with ellipticity 
constants and 5. Then, for all X G C, L + XI is a continuous operator from W l,2 (T) 
to W~ 1,2 (T). Moreover, when X = 1, L + I is an isomorphism from W 1,2 (T) onto W~ 1,2 {T) 
and the norm of its inverse is bounded by a constant only depending on 5. 

Proof: That L + XI is bounded from W^^fT) to W~ 1,2 (T) was already seen in Remark 15.21 
above. The Lax-Milgram theorem, applied with the sesquilinear form 

B(u,v) := 2^ c x, y du(x, y)dv(x, y)fJ, xy + 2j u(x)v(x)m(x) 

which is clearly continuous and coercive on VP 1,2 (r) thanks to (E32D, yields the invertibility 
of L + / and the bound on the norm of the inverse. □ 

Relying on Lemma [4.41 again and arguing as in Proposition I4.5[ we prove: 

Proposition 5.6 Assume that V satisfies ( fDj) and (P2) and L is an elliptic operator as 
above. Then there exists e > only depending on the geometric constants of F and the 
ellipticity constants of L such that L + I is invertible from W l,p to W~ 1,p forp G (2 — e, 2 + e). 

5.4 Proof of the estimates on L 

We prove here Theorems 15.31 and 15.41 



5.5 Proof of Theorem 15.31 

Let / G £ 2 (T). Since L + I is an isomorphism from H /1 ' 2 (r) onto VF -1 ' 2 (T), there exists a 
unique u G W 1,2 (T) such that (I5.33j) is satisfied. Moreover, we have seen in Proposition 15.61 
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that, for some e > L+I is an isomorphism from W 1,P (T) onto W 1 ' P (T) for p G (2— e, 2+e). 
Let p G (2, 2 + e). Propositions 13.13 1 and 13.141 applied with q = a = 2 yield 

L 2 (r) w- 1 ' p (r) 

and 

^(r) m> c"(r). 

Therefore (L + J)" 1 maps L 2 (r) into C»(r) with 77 = 1 — |. 

With an analogous argument, L + I can be replaced by L + XI for A G with |A| = 1: 

IKL + Ar^lU^CH^D (5.34) 

and 

KL + A)- 1 ^ < C||n|| L2(r) . (5.35) 
We now claim that there is a constant C > such that, for all A G S M \ {0}, 

IKL + Ar^iu^qAi^iMi^r) 

and 

+ X)-^ <c\\\Z-$\\u\\vw 

We treat the case where A > using a scaling argument. A similar proof can be done for a 
general A G £ M . 

The scaling argument is as follows. For all a > 0, let T Q be the set T, endowed with the 
weights ah xy and the measures o?[i xy . Note that, if V Q stands for the gradient in r a , one 
has V a f — -V/. As a consequence, it is easily checked that, endowed with the distance d a 
defined by the weights ah xy and the measure a 2 m, the graph T a satisfies the assumptions 
(D), (P 2 ) and (L 2 ) with the same constants as T. 

Consider now the operator L a on T a given by the coefficients c xy . A straightforward com- 
putation yields that L a has the same ellipticity constants on T a as L on T. 
Let / G L 2 (T) and A > be given. Define g := jf. Applying the conclusion of Theorem 15.31 
with and the operator lyr, one obtains that there exists a unique function u G W 1,2 (T^) 
such that 

( L Vx + I ) u = 9 

on T^, and 

C 

IMW^) < C WgWvp^) = ^ ll/llw(r) • 
Therefore, (L + A/) it = / on T. Moreover, 

C 



\ u \\L°°(r) — II-' Hl-'(L'i 



and 

77-1 

"U U,r™ < CX 2 



\Cv(T) — ^ A IM llL 2 (r) 



which ends the proof of Theorem 
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5.6 Proof of Theorem 15.41 

The proof of Theorem 15.41 is analogous to that of Proposition 28, p. 59 in [1] in the Euclidean 
case, and we only indicate the main differences. 
Let us first recall that, for all t > 0, 

e~ tL = — I e zX (L + A) _1 dA (5.36) 
27U J 1 

where 7 is made of two rays y± := {re ±t9 ; r > 7} and of the arc 70 := {\e ta ; \a\ < 6} where 
g]|,7t[, and 7 is described counterclockwise. 
Combining (I5.36P with Theorem \5.'S\ we obtain 

He^itHL^D < Ct^^\\u\\ L 2^, t > 

and 

\e~ tL u\ dv < Ct~^ + ^\\u\\ 2 , t>0. 
We finish the proof as in Chapter 5, section 5.2, see also [T], Chapter 1. □ 
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